Quantum Mechanics helps in searching for a needle in a haystack 
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Quantum mechanics can speed up a range of search applications over unsorted data. For example imagine a phone directory con- 
taining N names arranged in completely random order. To find someone's phone number with a probability of 50% , any classical 

algorithm (whether deterministic or probabilistic) will need to access the database a minimum of 0.5N times. Quantum mechan- 
ical systems can be in a superposition of states and simultaneously examine múltiple names. By properly adjusting the phases of 
various operations, successful computations reinforce each other while others interfere randomly. As a result, the desired phone 

number can be obtained in only 0( JN) accesses to the database. 

1. Introduction In 1994 Shor discovered a quantum 1.1 Quantum Mechanical Algorithms In a quantum 



mechanical algorithm for factorization that was expo- 
nentially faster than any known classical algorithm [1]. 
This paper presents a quantum mechanical algorithm for 
search that is only polynomially faster than any classical 
algorithm; however, it does not depend for its impact on 
the unproven difüculty of the factorization problem. The 
search problem is this: there is an unsorted database 
containing N items out of which just one item satisfies a 
given condition - that one item has to be retrieved. Once 
an item is examined, it is possible to teli whether or not 
it satisfies the condition in one step. However, there does 
not exist any sorting on the database that would aid its 
selection. The most efficient classical algorithm for this 
is to examine the items in the database one by one. If an 
item satisfies the required condition stop; if it does not, 
keep track of this item so that it is not examined again. It 
is easily seen that this algorithm will need to examine an 
average of 0.5N items before finding the desired item. 

It is possible for quantum mechanical systems to 
make interaction-free measurements by using the dual- 
ity properties of photons [2]. In these the presence (or 
absence) of an object can be deduced by allowing a 
small probability of a photon interacting with the object. 
Therefore, most probably the photon will not interact, 
however, just allowing a small probability of interaction 
is enough to make the measurement. Thus in the search 
problem also, it might be possible to find the object 
without examining all of the objects, but just by allow- 
ing a certain probability of examining the desired object. 

Indeed, this paper shows that by using the same 
amount of hardware as in the classical case, but by hav- 
ing the input and output in superpositions of states, we 

can find an object in 0(jN) quantum mechanical steps 
instead of O(N) classical steps. Each quantum mechan- 
ical step consists of an elementary unitary operation 
(discussed in the following paragraph). 



computer, the lògic circuitry and time steps are essen- 
tially classical, only the memory bits that hold the vari- 
ables are in quantum superpositions (see [1] & [3] for a 
more detailed introduction to quantum computers). 
Quantum mechanical operations that can be carried out 
in a controlled way are unitary operations that act on a 
small number of bits in each step. The quantum search 
algorithm of this paper is a sequence of such unitary 
operations on a pure state, followed by a measurement 
operation. The three elementary unitary operations 
needed are the following. First is the creation of a super- 
position in which the amplitude of the system being in 
any of the N bàsic states of the system is equal; second 
is the Walsh-Hadamard transformation operation and 
third the selective rotation of the phases of states. 

A bàsic operation in quantum computing is the 
operation M performed on a single bit that is represented 



i.e. a bit in 



by the following matrix: M = — — 

the state is transformed into a superposition in the two 
states: \-^=, —=] . Similarly a bit in the state 1 is trans- 



formed into 



— — I, i.e. the magnitude of the 

J2 JïJ 



amplitude in each state is — but the phase of the 

amplitude in the state 1 is inverted. The phase does not 
have an analog in classical probabilistic algorithms. It 
comes about in quantum mechanics since the ampli- 
tudes are in general complex. In a system in which the 

states are described by n bits (it has N = 2 n possible 
states) we can perform the transformation M on each bit 
independently in sequence thus changing the state of the 
system. The state transition matrix representing this 

operation will be of dimension 2 n X 2". In case the ini- 
tial configuration was the connguration with all n bits in 
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the first state, the resultant configuration will have an 

1 

identical amplitude of 2 in each of the 2" states. This 
is a way of creating a superposition with the same 

amplitude in all 2 n states. 

Next consider the case when the starting state is 

another one of the 2" states, i.e. a state described by an n 
bit binary string with some Os and some ls. The result of 
performing the transformation M on each bit will be a 
superposition of states described by all possible n bit 
binary strings with amplitude of each state having a 

n 

magnitude equal to 2 and sign either + or -. To 
deduce the sign, observe that from the definition of the 



matrix M, i.e. M = 4= 

72 



1 1 
1 -1 



, the phase of the result- 



ing configuration is changed when a bit that was previ- 
ously a 1 remains a 1 after the transformation is 
performed. Hence if x be the n-bit binary string describ- 
ing the starting state and y the n-bit binary string 
describing the resulting string, the sign of the amplitude 
of y is determined by the parity of the bitwise dot prod- 

uct of x and y, i.e. (-1)* y . This transformation is 
referred to as the Walsh-Hadamard transformation [4]. 
This operation (or a closely related operation called the 
Fourier Transformation) is one of the things that makes 
quantum mechanical algorithms more powerful than 
classical algorithms and forms the basis for most signifi- 
cant quantum mechanical algorithms. 

The third transformation that we will need is the 
selective rotation of the phase of the amplitude in certain 
states. The transformation describing this for a 2 state 



system is of the form: 
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, where j = J-ï and 



(^p §2 ^ arbitrary real numbers. Note that, unlike the 

Walsh-Hadamard transformation and other state transi- 
tion matrices, the probability in each state stays the 
same since the square of the absolute value of the ampli- 
tude in each state stays the same. 

2. The Abstracted Problem Let a system have 

N = 2 n states which are labelled S 7 ,S 2 ,-% These 2" 
states are represented as n bit strings. Let there be a 
unique state, say S v , that satisfies the condition C(S V ) = 
1, whereas for all other states S, C(S) = (assume that 
for any state S, the condition C(S) can be evaluated in 



unit time). The problem is to identify the state S v . 

This could represent a database search problem 
where the function C(S) is based on the contents of 
memory location corresponding to state S (as discussed 
in the abstract). Alternatively it could represent a prob- 
lem where the function C(S) was being evaluated by the 
computer. Various important computer science problems 
can be represented in this form [3] [5] [9]. 

3. Algorithm 

Steps (i) & (ii) are a sequence of elementary unitary 
operations of the type discussed in section 1.1. Step (iii) 
is the final measurement by an external system. 



i.e. there is the same ampli- 



(i) Initialize the system to the superposition: 
J_ J 1_ _M 

-Jn' Jn' Jn'" JnJ 

tude to be in each of the N states. This superposition 
can be obtained in 0(\ogN) steps, as discussed in 
section 1.1. 

(ii) Repeat the following unitary operations 0{Jn) 
times (the precise number of repetitions is impor- 
tant as discussed in [5]): 

(a) Let the system be in any state S: 

In case C(S) = 1 , rotate the 
phase by 71 radians; 
In case C(S) = , leave the 
system unaltered. 

(b) Apply the diffusion transform D which is 

defined by the matrix D as follows: 

D. = £ if i * j & D H = - 1 + J 

(D can be implemented as a product of 3 ele- 
mentary matrices as discussed in section 5). 

(iii) Measure the resulting state. This will be the state S y 
(i.e. the desired state that satisfies the condition 
C(S V ) = 1 ) with a probability of at least 0.5 . 

Note that step (ii) (a) is a phase rotation transformation 
of the type discussed in the last paragraph of section 1.1. 
In an implementation it would involve a portion of the 
quantum system sensing the state and then deciding 
whether or not to rotate the phase. It would do it in a 
way so that no trace of the state of the system be left 
after this operation so as to ensure that paths leading to 
the same final state were indistinguishable and could 
interfere. Reference [5] gives a way of doing this with a 
single quantum query. Note that this does not involve a 
classical measurement. 
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4. Convergence The loop in step (ii) above, is the 
heart of the algorithm. Each iteration of this loop 

increases the amplitude in the desired state by 0\ -|= | , 

as a resultin O(jN) repetitions of the loop, the ampli- 
tude and hence the probability in the desired state reach 
(9(1) . In order to see that the amplitude increases by 

0^-j=^j in each repetition, we nrst show that the diffu- 

sion transform, D, can be interpreted as an inversion 
about average operation. Just a simple inversion is a 
phase rotation operation and by the discussion in the last 
paragraph of section 1.1, is unitary. In the following dis- 
cussion we show that the inversion about average opera- 
tion (defined more precisely below) is also a unitary 
operation and is equivalent to the diffusion transform D 
as used in step (ii)(a) of the algorithm. 

Let oc denote the average amplitude over all states, 

■th 



a vector each of whose components is equal to the aver- 
age of all components. 

2 

Using the fact that P = P , it follows immediately 

2 

from the representation D = - 1 + 2P that D = I 
and hence D is unitary. 

In order to see that D is the inversion about aver- 
age, consider what happens when D acts on an arbitrary 
vector v . Expressing D as - / + 2P , it follows that: 
Dv = (-I + 2P)v = -v + 2Pv. By the discussion 
above, each component of the vector Pv is A where A is 
the average of all components of the vector v . Therefore 
the / th component of the vector Dv is given by 
(-v i + 2A) which can be written as (A + (A-v ; )) 

which is precisely the inversion about average. 

Next consider the situation in figure 2, when this 
operation is applied to a vector with each of the compo- 

C 



i.e. if OC, be the amplitude in the i th state, then the aver- nents, except one, having an amplitude equal to 



4n 



age is L y a . . As a result of the operation D, the 
N i—i 1 



i = 1 



amplitude in each state increases (decreases) so that 
after this operation it is as much below (above) a, as it 
was above (below) oc before the operation. 
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Figure 1. Inversion about average operation. 

The diffusion transform, D , is defined as follows: 
(4.0) D u = l,ifi*j&D u = -l + l. 

Observe that D can be represented in the form 
D = - 1 + 2P where / is the identity matrix and P is a 

projection matrix with P ^ = for all i, j . The follow- 
ing two properties of P are easily verified: nrst, that 

2 

P = P & second, that P acting on any vector v gives 



where C lies between - & 1 ; the one component that is 



different has an amplitude of -Jl - C 



Average 
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Figure 2. Inversion about average operation is 
applied to a superposition where all but one of the 
components are initially identical and of magnitude 

1 



O 



The average A of all components is approximately equal 
C 

to — . Since each of the (N - l) components is 

approximately equal to the average, they do not change 
significantly as a result of the inversion about average. 
The one component that was negative, now becomes 
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positive and its magnitude increases by 



2C 

Jn' 



In the loop of step (ii) of section 3, first the ampli- 
tude in a selected state is inverted (this is a phase rota- 
tion and hence a vàlid quantum mechanical operation as 
discussed in the last paragraph of section 1.1). Then the 
inversion about average operation is carried out. This 
increases the amplitude in the selected state in each iter- 
2C 

ation by — . Therefore as long as the magnitude of the 

Jn 



amplitude in the single state, i. e. J 
1 



1 - C , is less than 
1 



,_ , the increase in its magnitude is greater than 

72 J2N 



It immediately follows that there exists a number M less 

than «Jn , such that in M repetitions of the loop in step 
(ii), the magnitude of the amplitude in the desired state 

will exceed . Therefore if the state of the system is 

now measured, it will be in the desired state with a prob- 
ability greater than 0.5 . 

5. Implementa tion As mentioned in section 1.1, 
quantum mechanical operations which can be imple- 
mented in terms of elementary unitary operations, are 
local transition matrices, i.e. matrices in which only a 
constant number of elements in each column are non- 
zero. The diffusion transform D is defined in step (ii)(b) 



of the algorithm as: D 



D as presented above, is not a local transition 
matrix since there are transitions from each state to all N 
states. Using the Walsh-Hadamard transformation 
matrix (section 1.1), D can be implemented as a product 
of three local unitary transformations as D = WRW , 
where R the phase rotation matrix & W the Walsh-Had- 
amard Transform Matrix are defined as follows: 

R = if i * / ; /?,, = 1 if i = ; R„ = -1 if i * . 

ij j n n 

W í; = 2~ n/2 (-l)' '■> ' , i is the binary representation of 
i , and ;' • j denotes the bitwise dot product of the two n 

bit strings i and j . 

Each of W & R is a local transition matrix. R as 
defined above is a phase rotation matrix and is clearly 
local. W, when implemented as in section 1.1, is a local 
transition matrix on each bit. 

We evaluate WRW and show that it is indeed 
equal to D. R can be written as R = R l + R 2 where 



R l = -I , I is the identity matrix and R 2 00 = 2, 
R 2 = if i * 0, j * . By observing that MM = I 
where M is the matrix defined in section 1.1, it is easily 
proved that WW=I and hence D l = WR X W = -I. We 
next evaluate D 2 = WR 2 W. By Standard matrix multipli- 

cation: D 2 ad = ^ W ab R 2,bc W cd- Usin S the defini " 

bc 



tion of R 2 and the fact N = 2 , it follows that 
D, 



2, ad 



™ a0 w 0d 



A(_l)a·Ò + Ò·à = J Thus 
2 n N 



all elements of the matrix D 2 equal — , the sum of the 

two matrices D 1 and D 2 gives D. 

The quantum search algorithm of this paper is 
likely to be simpler to implement as compared to many 
other known quantum mechanical algorithms. This is 
because the only operations required are the Walsh-Had- 
amard transform & the conditional phase shift opera- 
tion, both of which are relatively easy as compared to 
operations required for other quantum mechanical algo- 
rithms [6]. Also, quantum mechanical algorithms based 
on the Walsh-Hadamard transform (e.g. the search algo- 
rithm of this paper, [4], [7], [8]) are likely to be much 
simpler to implement than those based on the "large 
scale Fourier transform" [1], [6]. 
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